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ABSTRACT 

"A  is  presented  with  the  property  that  the 

path  follows  the  hoaotopy  path  tracked  by  the  usual  vector  labeling.  Ibis 
produces  an  algorithm  which  deteralnes  the  path  without  pivoting,  on  a linear 
system  and  without  the  extra  dlmanslon  introduced  in  the  '”saadvlch*approach. 
The  results  should  therefore  lead  to  computational  savings. 


and  RSF  Gnat  Bo.  08-^2010  . 
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1.  IHTR0DUCTI05 


The  theory  of  aimplicial  approximation  has  been  applied  by  numerous 
authors  to  problems  involving  the  computation  of  fixed  points  and  the  equi- 
valent problem  of  finding  solutions  to  a system  of  nonlinear  equations  [l] , 

[2],  [3],  [k],  [5].  [6],  [T].  [8],  [91.  [10],  [11],  [12],  [13],  [Ik],  [15]. 

[16],  [17].  In  moist  of  the  previous  works  in  this  field  two  apparently  differ- 
ent approaches  can  be  reco^ized.  These  two  approaches  pivot  upon  the  same 
triangulation  (of,  for  example,  Rn  * [0,  l] ) but  they  differ  in  that  one 
uses  a "scalar  labeling"  and  the  other  a "vector  labeling."  The  scalar  label- 
ing was  initially  presented  by  Scarf  [15]  for  the  problem  of  computing  fixed 
points  on  the  unit  simplex,  and,  by  simple  extension,  for  a asp  f:  C C, 

C compact  in  8s . Scalar  labelings  were  then  presented  by  Fisher  and  Gould  [5a] 
for  the  nonlinear  complementarity  problem  and  by  Fisher,  Gould  and  Tolle  [6],  [7], 
and  Gould  and  Tolle  [12]  for  the  general  problem  of  solving  f(x)  ■ 0,  f:  R11  *►  8°. 
The  so-called  "vector  labelings”  were  developed  most  notably  by  Eaves  [2],  [3]. 

[k],  and  Merrill  [13],  both  for  the  fixed  point  problem  and  the  problem  of  solv- 
ing f(x)  ■ 0 on  bounded  and  unbounded  regions.  Theorems  regarding  the  existence 
of  honotop y paths  were  shown  in  Garcia  [9].  and  Charaes,  Garcia  and  Leake  [1]. 

Merrill  [13]  end  Eaves  [5]  independently  solved  the  restart  problem  for  vector 
labelings.  The  problem  was  solved  by  Fisher,  Gould,  and  Tolle  [7]  for  scalar 
labelings. 

In  this  paper  we  shall  for  convenience  present  the  entire  discussion  in 
the  framework  of  seeking  a solution,  say  x*,  to  the  system  f(x)  ■ 0.  Although 
the  scalar  and  vector  approaches  appear  to  be  different,  it  is  known  tbat  in  a purely 
formal  way  any  scalar  labeling  t can  be  encompassed  within  the  vector  fiamssnili 
as  follows.  Let  v*  be  any  vertex  in  the  trlangulatlon.  Associated  with  this 

!• 


Shoogh  this  transition  la  partly  formal  aad  not  of  computational  interest  It 
iraqpts  the  eccnreree  question  as  to  whether  or  not  rae tor  can  ho 

txpraasod  is  the  scalar  fr  answer k.  For  esanple,  the  "sandwich"  woctor  procedure 
of  Morrill  [13]  assigns  a label  A*  to  each  vertex  v*  of  H®  x [o,  l]  as  fallova. 

If  r1  c H®  x (0) 
if  v*  t H®  x {1} 

OilbC  these  labels,  the  vector  approach  to  solving  f(x)  ■ 0 pivots  en  the 
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trl angulation  through  successive  basic  feasible  solutionr  of  the  syst< 


[“I  • [:].  “ 


eRn) 


vhere  IQ  indices  vertices  in  R°  * {0}  and  Ix  indices  vertices  in  R11  x (l). 
That  is.  the  algorithm  proceeds  via  complementary  pivoting  to  generate  a se- 


quence of  slmpllces  FQ,  SQt  F^  8^  Fj,  Sj,  ....  vhere  F^  and  FJ+1  are 
facets  (n-slapllces)  of  the  a+l-slmplex  8^.  The  move  from  Fj  to  8^  is 
made  by  Joining  a single  vertex  to  the  n+1  vertices  of  F . The  rule  for 
selecting  the  vertex  to  bring  in  is  given  by  the  trlanguletloa  formula.  The 


move  from  8^  to  Vl  is  made  by  dropping  a vertex  from  Sj.  The  drop-out 
rule  is  provided  via  the  linear  system  by  pivoting  from  the  basic  feasible 
solution  associated  vith  Fj  to  the  one  associated  with  Fj+^.  It  is  known 
(Charaes.  Garcia  and  Lemke  [ l])  that  when  the  slapllces  of  the  trlangulation 
are  suitably  small  the  path  generated  by  the  above  vector  labeling  is  (in  the 
limit)  the  bomotopy  path  R(x,  t)  • 0 vhere  H(x,  t)  5 tf(x)  ♦ (1  - t)x» 
t c (0.  1].  x « R° . Without  loss  of  generality  ve  have  chosen  the  initial  point 
to  be  OCR0.  Thus  1(0.  0)  ■ 0 and  if  the  algorithm  converges  to-  (x*.  1) 
me  have  H(x*.  1)  * 0 ■ f(x*). 

The  principal  result  of  this  paper  is  that  a scalar  labeling  sad  a modified 
pivoting  procedure  have  been  discovered  for  generating  precisely  the  same  homotopy 
path.  In  this  sense  the  vector  labeling  can  be  cast  into  the  scalar  frsmevork.  The 
result  is  of  computational  significance,  for  it  eliminates  both  the  extra  "sandwich" 

dimension  and  the  need  to  pivot  on  the  linear  system. 

2.  IKS  ALGORITHM  AHD  THE  LABELIHG 

*R*t  fs  R°  ♦ R°  is  continuous.  f^O)  j*  0,  i ■ 1,...,  n,  f(x*)  * 0, 


and  x*  is  interior  to 


orthant.  Each  orthant  of  Rn  is  triangulated  (in 
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general,  orthant a tn  determined  by  a starting  point  ut  rather  the 
origin).  There  is  no  need  to  add  an  extra  dimension  by  triangulating  R11  * [0,  l]. 
The  algorithm  begins  in  some  orthant  denoted  at  the  origin  with  an  n— simplex 
(n  ♦ 1 vertices)  containing  the  labels  1 through  n + 1.  A simplex  with  labels 
1 through  n + 1 is  tensed  an  (n  + 1) -complete  (or  a completely  labeled)  simplex. 
Complementary  pirating  with  scalar  labels  produces  a path  in  EJ  of  n-aiapllcea 
with  labels  1 through  a.  A simplex  with  labels  1 through  n is  termed  an 
n-complete  (or  a (a  ♦ 1) -almost  complete  simplex) . We  note  that  although  an 
(a  + 1) -complete  simplex  must  be  an  n-s implex,  an  n-complete  simplex  may  be  an 
n-eimplex  or  an  (n  - 1) -simplex.  The  slmpllcial  path  which  we  generate  is  des- 
cribed by  the  following  possibilities. 

(Pi)  The  path  consists  of  an  unbounded  sequence  of  n-complete 
a-slJVliceo  la  / 

(P2)  The  path  terminates  at  an  (n  + l) -complete  simplex  which 
is  near  x*  (Interior  to  X^). 


(P3)  The  path  terminates  at  an  n— complete  (n  - l) -simplex  in  the 
boundary  of  X^. 

In  eaee  (F3)  a nodi  fled  pirating  technique  will  generate  a sequence  of  boundary 
(n  - l)-slaplices  la  such  a way  that  either 

(Ft)  (a)  the  boundary  path  Is  unbounded,  or 
(Ft)  (b)  the  path  reverts  back  into  the  orthant  E1  or  it  breaks 
throng  the  boundary  into  a new  orthant  X^. 

If  (P4)  (b)  occurs,  than  one  of  the  poesibilitlee  (Pi),  (P2),  or  (P3)  will 
again  occur,  with  X*  replaced  by  X^. 

It  will  be  shown  that  when  the  triangulation  grid  is  suitably  msall  the  above 
described  path  follows  the  some  hanotopy  path  which  is  followed  by  the  vector  labeling. 


* 
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A new  scalar  labeling  is  Introduced  for  our  purpose.  At  each  step 
of  the  algorithm  the  current  simplex  either  intersects  the  interior  of  an 
orthaat  or  lies  in  a boundary  of  an  orthant.  He  shall  tern  the  "current 
orthant"  to  be  either  the. orthant  whose  interior  the  current  simplex  inter- 
sects or,  in  case  the  current  sinplex  is  in  a boundary,  the  current  orthant 
la  the  one  whose  interior  was  last  visited. 

Let  t1  denote  the  current  orthant  and  suppose  x e t? . It  • 0 
for  sobs  1,  define 

«•  if  t±(x)  >0  and  j±  > 0 for  all 

y in  or  it  t±(x)  < 0 and  y4  < 0 

for  all  y . in 

If  f^x)  >.  0 and  yt  < 0 for  all 

y in  t*  or  if  f^(x)  < 0 and  y^  ^ 0 

for  all  y In  ^ 

(depending  on  the  orthant  t* ) is  specified 

xifl(x)  i °»  i 

ft(x)  f (x) 

< “ ■■  all  J,  otherwise 

xi  J 

3.  m rarriAL  simple: 

We  now  identify  a starting  n-slaplex  which  is  a unique  (n  + l)-conplete 
sinplex  in  a neighborhood  of  thu  origin,  nie  space  8°  is  triangulated'  in  such 


/ 


Mx) 


< 


\ 

Bow,  for  any  x c t*  our  labeling 
as  follows: 


(1) 


Aj(x) 


a ♦ 1 if 
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ft  way  that  thft  Interior  of  any  a-sisplex  win  not  intersect  any  hyper- 
plaae  Xj«0  for  any  1.  Thus  each  n-sl^plex  is  entirely  in  an  orthant 
of  l“. 

As  stated  earUar  It  Is  assueed  that  f^O)  i 0,  i ■ 1,  ....  n.  Let  D denote 
the  n a a aatrlx  whose  diagonal  tarns  are  given  by  ■ -sgn  f^O)  and 
whose  off-diagonal  tens  4^  ( J 4 i)  are  zero.  Let  e1  be  the  ith  unit 
vector  in  K11  and  a be  the  n^vector  of  ones.  Finally,  let  c > 0 be  the 
grid  slse  and  triangulate  in  such  a way.  as  to  include  the  singles  (t°,  v1,  ...»  v^} 
where 

T°.0 

v*  • c D(e  - e1),  i • 1,  ....  a . 

Thus,  for  1*1,  ...,  a,  ■ 0 and,  for  J f*  1,  v^  ■ -e  sgn  f^CO).  This 
a-slnples  will  be  the  starting  singles,  denoted  o£  and  it  will  be  shown  that 
this  is  a unique  (n  ♦ l)-conplete  singles  in  son  neighborhood  of  the  origin. 

Let  H1  denote  the  orthant  containing  ojj.  Vote  that  if  s is  in  the  interior 
of  ^ than  < 0 for  all  4. 

L— a 1.  Let  1(0,  d)  be  ft  neighborhood  of  the  origin  for  which  f^(s)  f 0, 
aU  1,  for  every  s in  1(0,  d).  Than  for  c g 0 sufficiently  snail  the 
singles  oj  is  a unique  (n  ♦ D-coaglete  sisglew  contained  in  V(0,  d). 

Proof-  It  is  lenadiate  fron  the  labeling  rules  (l)  that  ij(v°>  • n ♦ 1. 

Let  c be  chosen  so  Mall  that  oj  la  ooetalned  in  1(0,  d).  lota  that  for 
any  s in  I1  fl  1(0,  d)  such  that  s^  i 0 It  anst  bo  the  ease  that  s^  and 
fj (a)  have  opposite  aiM«-  Pick  an  index  U (l,  ...,  a).  leoaU  • 0 end 
Wj  f*  0,  4 4 i.  Hence, 


WM 


i 


or  a- . Otherwise,  If  r lith*  T«r*  of  <T  labeled  a ♦ 1,  than,  frcai  the 
laballa*  mla  (l),  WjfjtT  ) > 0,  all  j.  If  acae  r±  * 0,  than,  slnca  f (▼) 
lx  also  nonzero,  It  follows  that  both  htn  tha  same  sign  which  contradict*  tha 
fact  that  for  anjr  x la  H1  0 1(0,  3)  It  mist  be  that  If  xt  i 0 than  x l and 
ft(x)  hawa  opposite  signs.  Slnca  t°  (tha  origin)  is  In  0®,  there  la  a 
coordinate  1 such  that  for  each  of  tha  other  wartlces  the  1th  coordinate  is 
nonxaro  and  soae  other  coordinate  la  aero  (aaa  Figure  1).  Hence,  If  o®  + oj 
tha  label  1 la  excluded  flroai  o®,  a contradiction.  f 
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k.  THE  PATH  H Ef1 

Starting  with  the  (n  ♦ l)-cooplete  « implex  o£  in  the  orthant  2J, 
th*  algorithm  generates  an  (a  ♦ 1) -almost  complementary  path,  i.e.,  a se- 
«me«  of  distinct  n-slmplices  {ojj,  oj,  crj,  ...}  whore 

(1) 

(11) 

(HI) 

Klthor  the  path  of  slapllcee  1«  unbounded  or  else  it  stop*  with 

(a)  either  an  (n  ♦ l) -complete  simplex  o£,  or 

(h)  an  n-complate  (n  - 1) -simplex  on  some  boundary  • 0 
• on  I*. 

*■»  SOW  show  that  in  the  limit  this  path  of  slmpllces  in  Z*  coincides  with 
the  hemotopy  path  H(x,  t)  • 0 where  H(x,  t)  3 tf(x)  ♦ (1  - t)x,  t e (0,  l], 
lots  that  H(x,  0)  » 0 if  and  only  if  x ■ 0.  and  if  t j»  0 then 
H(xt  t)  • 0 If  and  only  If  f(x)  » Xx  for  soma  X < 0.  Thus,  for  x + 0 the 
homotopy  path  Is  described  by  the  known  condition  [l] 

f(x)  ■ Xx,  some  X < 0 . 

StSEfiJL  ^*rt  *jj  be  a sequence  of  grid  sizes  such  that  ♦ 0 and  let 
C denote  the  connected  closed  set  which  Is  the  curre  of  the  shore 

generated  path  of  slmpllces  in  . Then,  for  each  x c C,  f(x)  • Xx  for 


each  oj  la  in  H1 

each  oj  Is  n-cosqplete,  for  J > 0 

^ is  an  n-coaplete  (n-l) -simplex,  all  J . 


X < 0. 
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Proof.  Suppose  x C C end  x^  # 0,  Xj  + 0.  Then  since  each  simplex  is 
n-coaplete , and  by  continuity,  it  must  be  that  t^x)/x^  ■ fj(x)/xj  <_  0. 

Thus  there  is  a X < 0 such  that  ffc(x)  » Xx^  for  each  xfc  * 0.  Suppose 
■ 0.  Then  in  each  neighborhood  of  x there  is  a point  y in  / with 
label  i and  hence  at  which  f1(y)/y1  <0.  In  each  neighborhood  of  x there 
is  also  a point  x c E*7  with  label  J + i.  Thus  £ tJLzVx^  Since 

ia  near  zero,  and  the  possible  values  of  Zj  can  be  bounded  away  from  zero, 
the  last  inequality  means  that  fl(z)/zJL  > 0.  Since  y and  z are  both  ia  E*7, 
it  aust  be  the  case  that  and  zt  have  the  same  sign.  Hence  f±(r)  and 
ft(*)  have  opposite  signs.  It  follows  that  f^x)  ■ 0 and  hence  x±  satisfies 
f^x)  • Xx^  for  any  X.  § 

Vote:  Hy  reasoning  analogous  to  that  in  the  proof  of  Theorem  1 it  can  be 

seen  that  any  point  In  a "small"  4 a ♦ 1 ) -complete  terminal  simplex  is  an  approxi- 
mate zero  of  f . 

5.  THE  ANALYSIS  WHET  THE  PATH  TERMINATES 
IN  THE  BOUNDARY  OP  EJ 

For  simplicity  of  exposition,  we  shall  make  the  nondegeneracy  assumption 
that  except  for  the  origin  the  homotopy  path  {x:  f(x)  ■ Xx,  X ^ 0}  intersects 
any  boundary  xt  • 0 only  at  points  x such  that  x^  i 0,  J i i. 

®°IW  aaawme  that  the  slmplicial  path  terminates  with  an  n-camplete  (n  - l)-slm- 
plex  in  the  boundary  of  E7 . Under  the  nondegeneracy  assumption  there  will  be  a 
E > 0 such  that  if  the  grid  c is  sufficiently  small  then  for  every  x in  the 
simplex  |xj  | £ C,  all  J t i.  A modified  pivoting  procedure  will  now  be  described 
for  use  on  x^  ■ 0 and'  it  will  be  shown  that  this  path  in  the  boundary  of  E7 
rams  Ins  arbitrarily  close  to  the  homotopy  path. 
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Suppose  0°  ^ is  an  n-conplete  simplex  in  that  boundary  of  for 
which  xt  ■ 0.  Let  a®-1  ■ {t*,  ...»  ▼“},  Hv*)  ■ J,  ■ 0,  all  J,  ' 
and  for  each  4 we  have,  by  nondegeneracy,  |v^|  > K all  m »*  i and  some  K > 0.  The 
(n  - 1) -simplex  an~1  is  a facet  of  exactly  two  n-simplices . One  of  these, 

A 

o* » it  in  . The  other,  d® , is  in  the  orthant  on  the  "opposite 
side"  of  the  boundary  xl  ■ 0.  Let  Ij  denote  the  labeling  function  on  the 
orthant  , where  the  rules  for  ij  are  also  given  by  (l) . Points  an  the 
boundary  x^  ■ 0 will  have  different  labels  depending  on  whether  the  labeling 
function  tj  or  is  employed.  This  is  because  for  such  points  f^(x)/x^ 
will  be  -*«•  on  one  side  of  ■ 0 and  -»  on  the  other  side,  depending  on 
the  sign  of  ^ it  x. 


l£5HS_2.*  Suppose  r is  in  the  boundary  x^  » 0,  where  lj.(v)  + n + 1 
(ii  tjlr)  - i i-  (T)  + i. 

(ii)  If  tj(w)  * 1 then  ij(v)  ■ smallest  J j*  i such  that 
*j(t)/Wj  <,  fk(T)/wk,  all  k i*  i. 


Proof.  Observe  that 


ij(T)  • i 


Bence,  (1)  holds,  and  (ii)  follows  immediately.  # 

Bow  let  o®"1  be  an  arbitrary  (n  - 1) -simplex  on  the  plain  x^  ■ 0 and 
let  a11  and  o®  be  the  n-slmpllces  in  X^  and  X^  respectively  such  that 

O®-1  - O®  n S®. 


(i)  O®"1  is  an  n-coaiplete  facet  of  a11  (with  respect  to  lj) 
iff  ijOJ®"1)  ■ U,  ....  n). 

(il)  <7U_X  is  so  n -complete  facet  of  a11  (with  respect  to  If) 

iff  4j(t)  i i for  a unique  v of  o0”1  and  Aj(ca“1  - {▼>)  » (l,  ...»  n)  - {i>. 

Proof.  Part  (1)  is  definitional.  Por  part  (ii) , a11'1  is  n -complete  (with 
respect  to  If  iff  ■ (1,  ...»  n}.  The  conclusion  follows  from  (i)  in 

Lane  2.  # 

Definition  1:  Let  a be  a simplex  on  » 0.  This  simplex  is  termed  jointly 

a-complcto  (jnc)  if 

(i)'  ij(r)  * i for  some  vertex  ▼ of  a. 

(ii)  lj(v)  ft  1 for  same  vertex  v of  a. 

(iii)  ij(o) U lj(o)  - {lf  ...»  n>. 

Vote  that  only  an  (n  - 2) -simplex  or  an  (n  - l) -simplex  on  x^  ■ 0 can  be  Jnc. 

Also*  by  Lemma  2,  tj  can  be  replaced  by  If  in  (i)  and  (ii>  of  the  above  definition. 

• • 

Lets  U.  Let  0°”^  be  a simplex  on  x^  • 0 such  that  0*-1  is  n-complete  with 
respect  to  either  lj  or  If.  Then  a11-1  bee  a Jnc  f«cet  o*  . 

Proof.  Suppose  or0”1  ■ {r\  ...»  v“).  Suppose  lJ(on“1)  • (it  ••••  b)  and 
let  tj(v^)  * J,  all  J.  Suppoee  tf(r^)  ■ k.  By  Learn  2,  k |*  1.  Than 

a0*1  - (v*>  satisfies  (!)•  (11) • and  (iii)  of  Definition  1 and  hence  is  Jnc.  Let 

<r  ■ a0"1  - (v^),  J ft  i,  k denote  any  other  facet  of  o°”1.  Them  J t ij(a)U  If  (a) 
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and  hanca  0®"1  - {v1}  la  not  Jnc.  Also,  a0"1  - {v1}  is  not  Jnc  because 
it  contains  no  vertex  r with  lj(v)  • i.  Hence  a11"1  - {vk}  is  the  only 
Jnc  facet  of  0®”1.  The  sane  argument  applies  if  i^c®-1)  ■ {l,  ....  n}.  # 


Lawns  5-  Let  cr11-1  be  a Jnc  simplex  on  » 0 but  not  n-coaplete  with 

respect  to  either  lj  or  ij.  Then  a®”1  has  exactly  two  Jnc  facets. 

Proof.  Let  0®”1  ■ {t1,  ...,  v®}.  Since  o®*1  is  Jnc  we  can  assume 

Al.  ^(r1)  - i,  ^(r1)  - k i i. 


A2.  for  soaw  s f*  1,  ij (▼*)  * s (s  could  equal  k) 
A3.  {Aj(v®),  ■ {j,  i},  for  J j*  i. 


Case  (i).  ij(v*)  ■ i.  Vote  that  v*  c (J®-1  - {r1)  and  r*  t a®"1  - {r1}. 

P»ts  (i)  and  (ii)  of  Definition  1 are  satisfied.  From  A3  plus  the  assunption 
that  lj(v*)  ■ i it  follows  that  0®"1  - {r1}  is  Jnc. 

Case  (ii).  t^v*)  - k i i.  There  is  a vertex  v e 0®"1  - {v1,  v*}  such 

that  tj(v)  ■ i.  Otherwise  l^o®"1)  • {l,  ...»  n}  which  is  not  possible  by 

the  hypothesis  of  the  la mm.  Since  ▼ e o®”1  - {r1}  and  v^  e O®”1  - {v*>, 
it  follows  aimla  that  tf®"1  - {v1}  is  Jnc. 

Bow  consider  the. facet  o®-1  - {v*}. 

Case  (i).  Ijtv^)  - i.  lota  that  v1  c 0®*1  - {v^}  and  v*  c O®"1  - {v*}. 
From  A3  and  the  fact  that  ijfr1)  ■ i it  follows  that  o®"1  - (v^)  is  Jnc. 

Case  (ii).  Ij(t^)  • k + 1.  There  is  a vertex  v c cr®”1  - {▼*,  vk}  such 

that  tj(v)  f*  i.  Otherwise  i^o®”1)  ■ {l,  ...t  n),  which  is  not  possible  by  the 
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hypothesis  of  the  Iona.  Since  v e o””1  - {rk>  and  r*  e a11-1  - {vk} , 
it  follows  that  0°  ^ - {▼*}  is  Jnc.  Hots  that  these  are  the  only  two 
Jnc  facets  of  a11”1.  § 


now  show  that  la  the  Unit  the  simplicial  path  in  the  boundary  of 
coincides  with  the  hoaotopy  path. 


Zet  x be  a point  of  convergence  of  a sequence  of  jnc  slapllces 
f(x)  • Xx  for  lose  X < 0. 


• 0. 


Ut  x be  a point  of  conrercence  of  a sequence  of  Jnc  slapllces 
2 on  xi  • 0.  hr  (1)  In  Definition  1, 

fjCx)/^  • on 

(11)  of  Definition  1, 

f^xl/x^  ^ 

oe  ft(x)  • 0.  Dr  (ill)  of  Definition  1, 

fj(x)/Xj  < f^tx)^  for  ere ry  k * 1,  J * i. 
ee  f(x)  • Xx  for  seas  X < 0.  $ 


oar  description  of  the  nodi  fled  pirating  in  the  boundary  of  E^. 
«rl-  the  first  n-cosvletd  facet  en  x^  ■ 0,  we  generate  a 
of  distinct  (n  - l)-siaplices  on  x^  ■ 0.  The  sequence  la  of  the  form 


t-S’1.  or1,  oj-1. ...» 


vtMT*  oj  1 n “j+^  nr.  Joe  for  oil  J , until  vo  obtain  a tamlnol  facet 
«rl  which  Is  again  an  n-ccnplete  facet  of  ^ or  / (ia  which  case  we 
either  pivot  "back"  Into  s or  "through"  the  boundary  Into  t*) . lots  that 
because  of  the  assumption  that  a solution  r*  Is  Interior  to  sosm  orthaat, 
for  a snail  grid  slse  we  cannot  encounter  a vertex  labeled  n ♦ 1.  Also,  by 
the  nondegeneracy  assumption,  we  cannot  intersect  the  boundary  of  the  pi 


x^  ■ 0 (for  a 


grid  size).  Voncycllng  Is  assured  as  cr^”1  and  «rl 
a jl-l  jn-1  J-l  w . 


have  only  one  Jnc  facet  whereas  oj  , o®  , ....  each  have  exactly  two 


such  facets. 


In  sunaary,  we  have  the  following  rules  for  pivoting  In  the  boundary  of 
on  which  x^  » 0: 

1.  first  (n  - l) -simplex,  o®"1,  Is  n-conplete.  Let  !j(v*)  ■ 1, 
^(v1)  ■ k.  Resort  the  vertex  s with  lj(x)  ■ k to  obtain  a®"2  ■ O®”1  - {t}, 
where  o®”2  1 a-  Jnc. 

2.  Given  any  o®“2  which  is  Jnc,  Join  a vertex  v by  the  trlangn- 
lation  rule  to  obtain  a new  o®"1  » c®~2U  {▼). 

a.  If  o®’1  Is  casuists  with  respect  to  Ij.,  return  Into  t* 
by  editing  an  (n  ♦ l)*^>  vertex  In  the  Interior  of  and  forming  a o®  In 

b.  If  o®-1  is  eoaplete  with  respect  to  ^(l^(v)  j*  1 for  a 
unique  v and  ti(a®“1  - {v})  ■ {l,  ....  n)  - {i})  go  through  x,  ■ 0 Into 
r by  adding  an  (a  ♦ 1)^  vertex  In  Interior  of  t1  end  fomlng  a 3®  in 

A 

** 

e.  Otherwise,  either:  case  (1)  lj(v)  • 1,  t~(v)  • p or 
case  (11)  lj(v)  ■ p#  1.  In  either  case,  raaove  the  vertex  t t*  v such  that 
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(a)  ijU)  * p,  or  (b)  Ijtz)  ■ i,  i~(x)  « p.  Since  on~2  io  Jnc 
either  (a)  or  (b)  vill  occur,  but  not  both. 


6.  GEOMETRIC  ILLUSTRATIONS 

Prior  to  thia  work  It  haa  been  quite  difficult  to  illuatrate  the 
courae  of  the  hcootopy  path.  Thia  atate  of  affaire  la  nov  considerably 
lnproved.  Aa  an  illuatrat ion , Figure  2 shovs  the  nev  scalar  labeling  for 
the  functions 

y)  • x - y - 10 
f2(x,  y)  * 2x  - y - 5 . 

Note  that  the  path  does  not  converge.  If  the  function  f vere  to  be 
replaced  by 

«<x)  • J^Oj-fU)  ■ 

then  a representation  slnilar  to  Figure  3 would' be  obtained  the  honotopy 
path  would  converge.  This  idea  of  using  the  Jacobian  transfornatlon  for  scalar 
label Inga  was  first  presented  by  Fisher,  Gould  and  Tolle  [ 7 ] and  was  later 
studied  by  Wolsey  [17]  end  Todd  [16].  The  present  work  underscores  the  impor- 
tance of  this  transformation  for  the  nev  scalar  presented  herein  mr,A 

hence  for  the  vector  'laK*Mng  as  well. 
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